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^ ■ Abstract: Models with essential non-perturbative QCD dynamics for the 

cn ■ 

^ . production at the LHC. Differences of up to an order of magnitude occurs, al- 

though the different models all describe the Tevatron data on high-p^ charmo- 
nium. An important point is here the treatment of higher order perturbative 
QCD effects. We also discuss the large rate of prompt J/ from these models 
background to CP violation studies. 

'f: : 1 Introduction 

> 

■ The important interplay of hard, perturbative QCD (pQCD) and soft, non-perturbative 
QCD effects has recently been clearly demonstrated in the production of heavy quarkonia 
in hadron collisions. The Tevatron data [||] for high-p_|_ J/ip, ip' and T is up to factors of 
50 above the pQCD prediction in the Colour Singlet Model (CSM) [0, where a colour 
singlet cc pair is produced at the parton level and gives a charmonium state with the 
same quantum numbers. This deficit can be explained by letting a fraction of the more 
abundant cc pairs in a colour octet state be transformed to a singlet state through some 
soft, non-perturbative QCD dynamics. Lacking a proper understanding of non-pQCD, 
this has been described in different models: the Colour Octet Model (COM) 0, the 
Colour Evaporation Model (GEM) [Q, the Soft Colour Interaction model (SCI) and 
the Generalized Area Law model (GAL) 0. All these models can be made to fit these 
Tevatron data. 

In this paper we study the extrapolation of these models to the LHC energy and 
examine the theoretical uncertainty in the charmonium production rate. Future LHC 
data may then discriminate between the models and provide an improved understanding 
of the non-pQCD mechanism they involve. 
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Prompt J / ip production is also of importance as a background to CP- violation studies 
based on B meson decays into J/ipX, with X being K^, etc. The B production cross 
section at the LHC is very high in comparison with e+e~ colliders, but its fraction of 
the total inelastic cross section is small (0.7%). The trigger for B physics has therefore 
to be very selective and is typically based on the leptons from a J/ip decay. Studies of 
such triggers have shown that prompt J/ip is indeed a source of background 0. The 
prompt J/ip, which emerge from the primary interaction vertex, can to a large extent 
be distinguished from J/tp from B decays, which are produced at a secondary vertex 
typically located a few hundred micrometers from the primary vertex. There is, however, 
a remaining prompt charmonium background that affects i?-physics analyses and the 
associated CP-violation studies, such as measurements of sin2(3. 

In Section 2 the models and their normalisation to the Tevatron data are discussed. 
The extrapolation to the LHC energy is made in Section 3, where also the observed 
differences are analysed and prompt J/ip as a background for CP- violation studies is 
considered. We end with some conclusions in Section 4. 



2 Models and tuning to Tevatron data 

In the Colour Octet Model the cross section is factorised in a short distance part, where 
a cc pair is produced in a well defined quantum number state (^'^+^Lj), and a long 
distance part, giving the probability that this state will convert non-perturbatively into 
a charmonium state. These probabilities are given by non-relativistic QCD (NRQCD) 
matrix elements, which are assumed to be universal and in practice are free parameters 
obtained from fits to experimental data. For high-p^ J/ '4^ production at the Tevatron, 
the main subprocesses are gg —>■ J/ipg and gq — J/ipq, which are next-to-leading order 
in the hard pQCD cc production process. The extraction of the NRQCD matrix elements 
from the Tevatron data has been performed in several steps where perturbative effects and 
intrinsic transverse momenta have successively been taken better into account explicitly 
instead of being absorbed into the fitted matrix elements 

The CEM, SCI and GAL models are based on a similar phenomenological approach, 
where soft colour interactions can change the colour state of a cc pair from an octet to a 
singlet. They employ the same hard pQCD processes to produce a cc pair regardless of its 
spin state. The leading order (LO) processes are gg —>■ cc (Fig. |I|a) and qq — > cc. Heavy 
quark production is, however, known to have large contributions from next-to-leading 
order (NLO) diagrams 0. Virtual corrections to the leading order processes together 
with soft and coUinear gluon emissions give an increase of the cross section, which can 
be approximately described as an overall i^'-factor multiplying the leading order cross- 
section. The LO processes and these NLO processes cannot produce a J/ip at high-pj_ 
since there is nothing to compensate its p± to give the essentially zero p± of the initial 
partons. Of importance for high-p^ J/ ^ production is instead NLO tree diagrams with a 
third hard parton that balances the p± of the cc pair. The most important contribution 
is given by the diagram in Fig. lb. Although this is an 0{as) correction in terms of a 
gluon sphtting g —>■ cc applied to the basic 2 — > 2 process of gg ^ gg, it is numerically 
large since gg — > gg has a much larger 0{a'^) cross section than the LO cc production 
processes. Matrix elements with non-zero quark masses are only available up to NLO. 
Still higher orders can be expected to be important at the Tevatron and LHC. The reason 



is that at these energies many gluons can be emitted and their virtuahty need not be 
very large in order to allow a split into a cc pair. The higher order processes can be 
approximately described by the parton shower approach available in the Pythia [|ry] 
Monte Carlo, where in all basic QCD 2 — 2 processes the incoming and outgoing partons 
may branch as described by the DGLAP equations 
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Figure 1: Illustration of cc production processes in (a) leading order (a^) and (b) next- 
to-leading order (af). 



In CEM the exchange of soft gluons is assumed to give a randomisation of the colour 
state. This implies a probability 1/9 that a cc pair is in a colour singlet state and produces 
charmonium if its mass is below the threshold for open charm production, rricc < 2m o- 
The fraction of charmonium giving a J/ip is given by an additional non-perturbative 
parameter pj/^ = 0.43 — 0.5 This model was recently implemented in Pythia such 
that higher order pQCD processes could be included in terms of parton showers and events 



be Monte Carlo simulated. It was found |12| that this model reproduced quite well xp 



and p± distributions, both in shape and normalisation, of J/ip produced in fixed target 
experiments. 

In the Soft Colour Interaction model |^ it is assumed that colour-anticolour, corre- 
sponding to non-perturbative gluons, can be exchanged between partons emerging from 
a hard scattering and hadron remnants. This can be viewed as the partons interacting 
softly with the colour medium, or colour background field, of the initial hadron as they 
propagate through it. This should be a natural part of the process in which 'bare' partons 
are 'dressed' into non-perturbative ones and the confining colour flux tube between them 
is formed. The hard parton level interactions are given by standard perturbative matrix 
elements and parton showers, which are not altered by softer non-perturbative effects. 
The unknown probability to exchange a soft gluon between parton pairs is given by a 
phenomenological parameter R, which is the only free parameter of the model. These 
colour exchanges lead to different topologies of the confining colour force fields (strings) 
and thereby to different hadronic final states after hadronisation. This model gives a novel 
explanation of rapidity gap events in deep inelastic scattering and in hard pp processes at 
the Tevatron |^, |13|, which are well reproduced with R = 0.5. Applying the same Monte 
Carlo implementation in Pythia (with the same i?- value), it was found that the Teva- 
tron data on high-p^ charmonium and bottomonium are also well reproduced |]14|. The 



increased production rate is here given by the possibility for a perturbatively produced 
QQ pair in a colour octet state to be transformed to a singlet state as a result of these 
soft colour interactions. The mapping of cc pairs, with mass below the threshold for open 



charm production, is here made based on spin statistics which avoids introducing further 
free parameters. This was also found to give a correct description of the different onium 
states observed at the Tevatron [Q. 



An ahernative to SCI is the later developed Generalised Area Law model [Q, where 
modified colour string-field topologies are obtained by interactions between the strings in 
the event instead of the partons. A generalisation of the area law suppression e~^^, with A 
the area swept out by the string in energy-momentum space, gives a dynamic probability 
R = Ro{l- e"''^^) for two string pieces to interact depending on the area difference 
AA resulting from the changed string topology. This favours making shorter strings and 
thereby favours quarkonium production. The parameter Rq and some hadronisation model 
parameters, e.g. b, are obtained from a fit to data from both deep inelastic scattering and 
e~^e~ annihilation 0. 

The comparison of the CEM, SCI and GAL models with the Tevatron data is shown 
in Fig. 1^. As can be seen, all models give a quite decent description of the data. Although 
the shape is not perfect in the tail of the distribution, it is quite acceptable given the 
simplicity of the models. The overall normalisation is correctly given by the models. For 
the CEM this is obtained by setting pj/^ = 0.43 and the charm quark mass to 1.5 GeV. 
As parameterisations of the parton density functions we have used CTEQ4L [|15|, but we 
have checked that the result does not change much if we use GRVHO or CTEQ2L {e.g. 
CTEQ2L reduces the normalisation by a few percent). The SCI and GAL models have 
not been tuned to these data, but the result is also sensitive to the charm quark mass 
(taken as default m,c = 1.35 GeV in Pythia 5.7). The parton densities are here kept the 
same as used in other applications of these models, namely CTEQ3L [|l6l for SCI and 



CTEQ4L [|T5| for GAL. The values of the essential parameters R and Rq in SCI and GAL 
are, as discussed, given by the rate of rapidity gap events in deep inelastic scattering. The 
SCI model is, however, quite stable such that R can be chosen in the range 0.2-0.5. One 
should note that an arbitrary K factor is not needed in any of the models, since higher 
order pQCD processes were included through the parton showers. 

In view of the Tevatron running with increased luminosity one can expect forthcoming 
J/ip data extending to higher p±, which would test and constrain the models further. We 
therefore include in Fig. § model predictions for p± up to 40 GeV, which show that the 
CEM and SCI models are quite close, whereas the GAL model gives a somewhat higher 
cross section in the high-p^ tail. New data could help in discriminating or improving the 
models resulting in reduced uncertainties in the model predictions. 

Based on the ability of these models to reproduce, in a reasonable way, the presently 
available Tevatron data, we now extrapolate them to the LHC energy. 



3 Extrapolation to the LHC energy 

Applying these models also for the larger energy at the LHC, i.e. y/s = 14 TeV, should 
be appropriate since they include a reasonable energy dependence. The production of 
the cc pair is given by hard pQCD processes with a well-defined energy dependence. The 
soft interactions that change colour octet states to singlet states have no explicit energy 
dependence, similar to the normal hadronisation process. The SCI and GAL models will, 
however, have an implicit energy dependence since they act on a parton state or string 
topology that depend on the collision energy. 
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Figure 2: Distribution in transverse momentum of prompt J/ip as observed by CDF and 
DO |T[] in pp interactions at the Tevatron and obtained in the CEM, SCI and GAL models. 



Results of the models, keeping all parameters fixed from the comparison with the 
Tevatron data, are shown Fig. ^. The p± distribution in Fig. |^a is quite similar for the 
three models, although they differ somewhat in the high-p^ tail and GAL has a somewhat 
less steep slope. These predictions should, however, not be taken as very precise in view of 
the simplicity of these models that attempt to describe unknown non-perturbative QCD 
phenomena. The overall normalisation, which e.g. is sensitive to the value of the charm 
quark mass, should not be considered to be better than within about a factor two. 

In order to consider what can realistically be measured, we show in Fig. §3 the dis- 
tribution in pseudorapidity of J/ip and the acceptance regions of the LHC experiments. 
Requiring the muons from the J/ ip decay to be within these regions we obtain the results 
in Fig. ^cd, which show lower effective cross sections with slightly more separated model 
curves. Thus, Fig. |^ shows a slightly different behaviour of the models as function of 
rapidity and transverse momentum, but a more detailed study of this is premature until 
some data are available to show that the models give a reasonable overall description of 
J/ip production at the LHC energy. A retuning of the models may then reduce these 
differences and it is therefore not clear to what extent these model differences can be 
exploited to discriminate among the models and disentangle details of the charmonium 
production mechanism. In any case, this would require data of quite high precision. 

Here we can make an interesting comparison with a previous calculation of the prompt 
J/ip production in the COM model As can be seen in Fig. |^c, the COM result is 
significantly lower than the other models. The COM result is based on using pQCD 
matrix elements in NLO. At high the dominant process is, as discussed above, the 
order 0{al) process gg ccg in Fig. |I|b. A shift of the momentum distribution of this 
cc pair has been made in order to account for the effect of higher order gluon emissions 
as estimated based on the parton shower. Since g ^ cc in the parton shower is here not 




Figure 3: Differential cross sections in transverse momentum and rapidity for J /ip in pp 
collisions at ^/s = 14 TeV based on the CEM, SCI and GAL models. In all cases only 
J/ip with p"^^ > 10 GeV and decaying into fi'^fi~, which in (c, d) are required to be 
within the indicated rapidity coverage of the Atlas and LHC-b detectors, respectively. 
For comparison, the COM results from [§] is included in (c). 

included, the applied higher order correction does not change the normalisation of the cc 
cross section, but only the shape of the momentum distribution which becomes somewhat 
softer due to gluon radiation. 

Our calculations with the CEM, SCI and GAL models use leading order matrix ele- 
ments and include the full parton shower evolution, including g —>■ cc in any branching. 
This gives an estimate of the cc production cross section including all higher orders. To 
find out to what extent this is the reason for the observed difference, we mimic the 0{al) 
matrix elements by including only those cc pairs coming from the first branching of the 
parton shower (cf. Fig. lb) giving the results shown in Fig. ^. Since we want to explicitly 
show the effect of omitting higher order cc production we have here not changed any pa- 



rameters in the models (which could be done to partly compensate for the loss of higher 
orders). As compared to the standard result, including cc from any branching, this 'first 
branching' approximation does indeed give a lower J/ip cross section. At the Tevatron 
energy the difference is not large and could at least partly be absorbed into a tuning of 
parameters such as and a^. Nevertheless, in the region of lower where the data are 
more precise, some preference for the all order result is indicated. At LHC, where more 
energy is available to build up a more extended parton shower, the difference is larger 
and becomes an order of magnitude at high p±. As can be seen in Fig. this reduction 
of the J /if) yield in the first branching approximation brings the CEM and SCI models 
closer to the COM result IQ. This COM result was based on the older parameterisation 
CTEQ2L of parton densities, but we have checked that this causes a reduction in the 
overall normalisation which is much smaller than omitting higher order g cc in the 
parton shower. Thus, we conclude that cc production in orders higher than 0{al) are 
important at the LHC energy. 




Figure 4: Distribution in transverse momentum of prompt J/ip at (a) the Tevatron and 
(b) LHC {\ri\ < 2.5 as for Atlas) for the colour octet model (COM) based on NLO 
0{al) matrix elements (from [Q) and for the colour evaporation model (CEM) including 
g cc to all orders in the parton shower and in the 'NLO' approximation where only 
cc pairs from the first branching have been included. The result of this 'first branching' 
approximation applied to the soft colour interaction model (SCI) is also shown in (b). 



As discussed in the Introduction, we also consider the prompt J/ip production as 
background for CP- violation studies based on B meson decays such as J/ipK^, 
-Bg J/'4"Pi 6tc. A detailed analysis has been made for Atlas fl^, taking into account 



trigger conditions, acceptance cuts, off-line selection criteria and reconstruction methods. 
This resulted in the signal-to-background ratio of 4.1 for the case of 5° — > J/ipcj). The 
estimated contamination from pp J/ipX on the level of 3% was here entirely due 
to prompt J/ip produced by the COM model implemented in the Atlas Monte Carlo 
package. Although other backgrounds were found to be more important, this signal-to- 



background ratio may be somewhat too optimistic in view of our finding that COM gives 
a lower prompt J/V' cross section than the CEM, SCI and GAL models. 

At this stage one cannot decide which of these models is most rehable and gives the best 
prediction for the prompt J / i]j production. The variation between them should therefore 
be taken as an estimate of the theoretical uncertainty. From Fig. 3c wc find cross sections 
that are almost an order of magnitude larger than in COM. The signal-to-background 
does not, however, decrease by the full amount of this factor since other backgrounds 
are also present. Nevertheless, this indicates that the prompt J/il) production could 
be an important background giving a lower signal-to-background ratio than previously 
estimated. A proper study of this is beyond the scope of this paper, but our results show 
that the prompt J/%Ij production is a background for CP- violation studies that should not 
be neglected. 

4 Conclusions 

Charmonium production provides an interesting testing ground for QCD, including both 
perturbative and non-perturbative effects. Several models (COM, CEM, SCI, GAL) have 
been developed to account for the possibility that a pcrturbatively produced cc pair in 
a colour octet state is transformed into a singlet state through soft interactions. This 
enhances the rate of charmonium production by more than an order of magnitude, as 
compared to conventional expectations based on the Colour Singlet Model. All these new 
models can account for the observed rate of high-p^ prompt J/%1) at the Tevatron, and are 
thereby normalised at this energy. Since the models tend to differ more at higher ■, we 
have extended our model calculations to give predictions for a region of higher p± which 
may be reached in the high luminosity runs at the Tevatron. Forthcoming high statistics 
data may then discriminate among the models or help in reducing the uncertainties in 
results of the models. 

Extrapolating the COM, CEM, SCI and GAL models to the LHC energy, we find 
significant differences in the predicted prompt J /il) cross sections; up to almost an order of 
magnitude. In particular, the COM result is lower than the others. Part of this difference 
is related to the fact that COM is based on NLO matrix elements, whereas the other 
models include still higher order cc production through the parton shower approximation 
in pQCD. 

Prompt J/'0 is also a background for studies of B meson decays into J/V', which are 
important for studies of CP-violation. An earlier estimate of the signal-to-background 
ratio based on COM gave the favourable result of 4.1. Using the larger rate of prompt 
J/tl) from the other models, will reduce the signal-to-background ratio. Prompt J/ if} 
production must therefore be better understood in order to control it as a background for 
CP-violation studies. 

Improving our understanding of the mechanism for prompt J/tl) production is also 
important in its own right. It involves an interplay of hard, perturbative and soft, non- 
perturbative QCD dynamics which is closely connected with the more general problem of 
understanding QCD. 
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